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LONG-TERM  GOALS 

Both  physical  oceanographic  processes  and  marine  geological  features  in  the  areas  of  continental  shelf 
and  shelfbreak  can  cause  horizontal  heterogeneity  in  medium  properties,  so  horizontal  refraction  and 
reflection  of  sound  can  occur  and  produce  significant  three-dimensional  (3-D)  sound  propagation 
effects.  The  long-term  goals  of  this  project  are  targeted  on  understanding  the  3-D  acoustic  effects,  and 
their  temporal  and  spatial  variability,  caused  by  the  environmental  factors  existing  commonly  in  the 
continental  shelf  and  shelfbreak  areas,  such  as  slopes,  submarine  canyons,  sub-bottom  layers,  surface 
waves,  internal  waves  and  shelfbreak  fronts. 

OBJECTIVES 

One  of  the  research  objectives  in  this  project  is  to  develop  efficient  and  accurate  numerical  models  for 
studying  underwater  sound  propagation  in  complex  ocean  environments.  The  ultimate  scientific 
objective  is  to  study  the  underlying  physics  of  the  3-D  sound  propagation  effects  caused  jointly  by 
physical  oceanographic  processes  and  geological  features.  To  achieve  this  goal,  individual 
environmental  factor  will  be  first  studied  and  then  considered  jointly  with  a  unified  ocean,  seabed  and 
acoustic  model.  Another  major  objective  is  to  develop  a  tangent  linear  model  to  predict  acoustic 
fluctuations  and  derive  sound  pressure  sensitivity  kernels  due  to  3-D  sound  speed  perturbation  in  the 
water  column.  The  sensitivity  analysis  will  consider  both  the  joint  ocean  and  seabed  effects. 

APPROACH 

The  technical  approaches  employed  in  the  3-D  sound  propagation  study  include  theoretical  analyses, 
numerical  computation,  real  data  analyses  and  field  work  experiments.  A  3-D  normal  mode 
method  has  been  used  to  study  canonical  environmental  models  of  shelfbreak  front  systems  [1]  and 
nonlinear  internal  wave  ducts  [2-3],  3-D  parabolic-equation  wave  propagation  models  with  improved 
split-step  marching  algorithms  [4-6]  are  used  to  study  sound  propagation  in  realistic  environments. 

The  parabolic-equation  (PE)  approximation  method,  first  introduced  by  Tappert  [1]  to  underwater 
sound  propagation  modeling,  has  long  been  recognized  as  one  of  the  most  efficient  and  effective 
numerical  methods  to  predict  sound  propagation  in  complex  environments.  The  advantage  of  this 
method  is  due  to  the  fact  that  it  converts  the  Helmholtz  wave  equation  of  elliptic  type  to  a  one-way 
wave  equation  of  parabolic  type.  The  conversion  allows  efficient  marching  solution  algorithms  for 
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solving  the  boundary  value  problem  posed  by  the  Helmholtz  equation.  This  can  reduce  significantly 
the  requirement  for  computational  resources,  especially  for  modeling  three-dimensional  (3-D)  sound 
propagation. 

A  higher  order  numerical  algorithm  has  been  proposed  to  split  the  square -root  Helmholtz  operator 
employed  by  the  PE  method  for  modeling  sound  propagation.  This  operator  splitting  method  deals 
with  multidimensional  cross  terms  to  yield  a  more  accurate  approximation,  and,  most  importantly,  it 
still  permits  efficient  3-D  PE  numerical  solvers,  such  as  the  Split-Step  Fourier  method  [5]  and  the 
Alternative  Direction  Implicit  (ADI)  Pade  method  [6].  The  higher  order  operator  splitting  algorithm  is 
shown  in  the  next  equation: 
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where  kref=  col cref  is  the  reference  wavenumber,  n  is  the  index  of  refraction  n  =  crejlc  =  k!krej ,  and 
denotes  the  Laplacian  operating  on  the  transverse  coordinates. 

Another  numerical  method  to  be  utilized  is  a  tangent  linear  solution  to  predict  acoustic  fluctuations  due 
to  3-D  sound  speed  perturbation  in  the  water  column.  This  model  has  been  developed  and  published 
[8],  and  it  is  briefed  here  along  with  an  application  to  determine  the  sensivitiy  kernel  of  sound  pressure 
variability  due  to  changes  in  medium  sound  speed.  Consider  the  one-way  parabolic  wave  equation, 

^cu{xry,z)  =  ikref|-l  + Jk-2fVi  +n2ay,z)j  u(x,y,z),  (2) 

where  u  =  p  exp (~ikTSlx),p  is  the  sound  pressure,  /cref  is  the  reference  wavenumber,  and  n  is  the  index  of 
refraction  with  respect  to  kre f.  Now,  denote  the  exact  PE  operator  as  the  following 


X  =  -1  +  Jkre2(V2±  +  n2  (x,y,z).  (3) 

and  let 

n2  (x,  y,  z)  =  y0  (x,  y,  z)  +  qq  {x,  y,  z),  (4) 


where  Yq  is  the  square  of  the  index  of  refraction  of  the  background  state,  and  hi  is  a  perturbation 
scaled  by  an  arbitrary  small  parameter  s.  Following  the  derivations  shown  in  [8],  one  can  find  the  next 
higher-order  tangent  linear  PE  solution,  which  unifies  two  other  less  accurate  tangent  linear  solutions 
proposed  by  Hursky  et  al.  [9]  and  Smith  [10]. 
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where  £°  1  +  skrefVj- +  Ya^x,y,z]  js  the  background  PE  operator. 

From  Eq.  (5),  we  can  deduce  the  sensitivity  of  sound  pressure  with  respect  to  n 2  to  be 
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and  with  the  chain  rule  we  can  also  get  =  -ikref  —  ^  (2  -  L0)p(x) .  This  sensitivity  kernel  is  in 

fact  the  local  kernel  at  the  perturbation  location.  To  obtain  the  kernel  at  the  receiver  location,  we  can 
multiply  the  local  kernel  with  the  Green  function: 
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WORK  COMPLETED 

The  tasks  completed  in  the  year  are  described  below. 

1.  3-D  sound  pressure  sensitivity  kernel 

From  the  higher-order  3-D  tangent  linear  PE  solution,  we  can  deduce  a  local  sound  pressure  sensitivity 
kernel  shown  in  Eq.  (6).  By  incorporating  the  Green’s  function,  we  can  then  extend  the  local 
sensitivity  kernel  to  the  receiver  position.  In  this  year,  an  adjoint  modeling  framework  has  also  been 
implemented  to  determine  the  dynamics  of  the  sensitivity  kernel  associated  with  the  physical 
oceanographic  field.  An  example  showing  the  sensitivity  kernel  of  sound  intensity  with  respect  to 
sound  speed  perturbation,  d\pfldc,  in  a  nonlinear  internal  wave  field  is  shown  in  Figure  1. 

2.  Joint  physical  oceanographic  and  acoustics  modeling 

A  framework  has  been  developed  for  integrating  numerical  models  of  ocean  dynamics  and  acoustics. 

A  data-assimilated  ocean  model  is  employed  to  drive  the  ocean  dynamics  in  this  joint  model.  An 
example  of  sound  propagation  on  a  shelfbreak  in  the  Mid- Atlantic  Bight  is  shown  in  Figure  2.  One  can 
see  that  the  original  down-slope  propagation  pattern  changes  to  water-borne  ducting  due  to  the  change 
of  water  column  stratification.  This  examples  demonstrates  the  use  of  this  integrated  numerical 
technical  for  studying  the  3-D  sound  propagation  effects  caused  jointly  by  physical  oceanographic 
processes  and  marine  geological  variability. 

3.  Engineering  test  in  the  New  England  Mud  Patch  area 

The  PI  also  joined  an  engineering  test  on  the  New  England  Mud  Patch,  south  of  Martha’s  Vineyard, 
MA  for  examing  propagation  conditions  of  broadband  acoustic  pulses  generated  by  a  combustive 
sound  source  deployed  by  ARL:UT.  One  of  the  hydrophone  arrays  deployed  in  this  test  was  provided 
by  the  PI  and  the  acoustics  group  at  WHOI.  Preliminary  acoustic  data  and  model  analyses  are  shown  in 
Figure  3.  Sound  propagation  models  in  two  different  sediment  types  of  environments  are  shown  in 
Figuers  3(a)  and  3(b).  The  sound  source  in  this  modeling  study  is  placed  close  to  the  seafloor.  In 
Figure  3(a),  the  sediment  layer  is  modeled  as  fine-granted  sediment  (mud)  with  compressional  wave 
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speed  less  than  the  sound  speed  in  the  bottom  of  water.  In  Figure  3(b),  the  sediment  layer  is  modeled 
as  sand  with  higher  compressional  wave  speed.  The  model  clearly  shows  that  the  acoustic  energy  can 
be  trapped  in  the  lower  speed  mud  layer.  A  preliminary  broadband  pulse  data  and  model  comparison  is 
shown  in  Figures  3(c)  and  3(d).  In  Figure  3(c),  the  spectrograms  of  received  pulses  on  the  hydrophone 
array  are  shown,  along  with  time  series  plots  in  the  bottom  panel.  In  Figure  3(d),  results  from  a 
broadband  sound  propagation  model  are  shown.  The  sediment  layer  in  the  model  is  made  to  be  mud. 
Despite  the  very  low  frequency  ground  waves,  this  preliminary  model  adequately  reproduces  the  pulse 
dispersion  over  a  broad  frequency  range. 

RESULTS 

The  major  results  of  this  project  accomplished  in  this  year  are  summaried  here,  along  with  a 
publication  list  provided  later.  First,  A  higher-order  tangent  linear  PE  solution  of  3D  sound 
propagation  has  been  derived,  and  it  unifies  other  tangent  linear  PE  solutions.  This  higher-order 
tangent  linear  solution  also  yields  a  sensitivity  kernel  of  sound  pressure  variations  due  to  changes  in 
medium  sound  speed.  An  adjoint  modeling  has  also  implemented  to  study  the  dynamics  of  the 
sensitivity  kernel  associated  with  the  physical  oceanographic  field.  A  numerical  example  of  3-D  sound 
pressure  sensitivity  analysis  in  a  nonlinear  internal  wave  field  is  presented  (Figure  1). 

A  numerical  model  has  also  been  constructed  to  include  the  sub-bottom  layering  structure  of  the  New 
Jersey  shelf.  The  purpose  of  this  model  is  to  investigate  the  3-D  sound  propagation  effects  caused 
jointly  by  marine  geological  features  and  oceanographic  dynamics  in  shallow  water.  Aligning  with  this 
research  topic,  the  PI  joined  an  engineering  test  on  the  New  England  Mud  Patch  for  a  preliminary 
sound  propagation  study.  The  initial  data  analysis  and  modeling  effort  is  presented  in  Figure  3.  In 
addition  to  that,  ambient  noise  analysis  has  been  initiated.  Figure  4  shows  the  noise  spectrum  level 
during  container  ships  passing  by  the  acoustic  engineering  test  site.  The  nominal  closest-point-of- 
approach  (CPA)  distance  to  the  hydrophone  array  was  about  5  to  7  km.  The  hydrophone  deployment 
was  5.5  days,  and  the  noise  spectrum  within  a  2  hour  time  window  around  every  CPA  are  analyzed. 
Total  35  containers  were  detected  acoustically,  and  the  total  analyzed  time  for  the  ambient  noise  was 
65  hours  long.  In  Figure  4,  the  black  lines  denote  the  mean  and  one  standard  deviation  curves.  The 
color  curves  are  the  average  levels  of  heavy,  moderate  and  light  shipping  noise  in  deep  ocean  [11],  The 
analysis  suggests  that  when  the  big  container  ships  pass  by,  the  ambient  noise  can  reach  up  to  1 10  dB 
re  lp  Pa  at  frequencies  around  100  Hz,  which  is  30  dB  more  than  the  average  level  in  deep  ocean 
reported  by  Urick  [1 1]  for  heavy  shipping. 

IMPACT/APPLICATIONS 

The  potential  relevance  of  this  work  to  the  Navy  is  on  increasing  the  capability  of  sonar  systems  in 
shallow  water  areas.  The  contributions  of  the  effort  on  studying  3-D  sound  propagation  effects  will 
also  be  on  assessing  the  environment-induced  acoustic  impacts. 

RELATED  PROJECTS 

Experimental  data  used  in  this  project  were  collected  from  the  ONR  SW06  and  QPE  projects  and  one 
of  the  the  reconnaissance  cruises  on  the  New  England  Mud  Patch.  In  addition,  strong  collaborations 
have  been  initiated  and  continued  between  the  Pi’s  of  an  ONR  MURI  project  on  integrated  ocean 
dynamics  and  acoustics. 
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Figure  1:  Sensitivity  analysis  of  an  internal  wave  acoustic  duct. 

[(a)  internal  wave  displacement  (the  source  is  underneath  the  first  wave  crest),  (b)  the  sensitivity 
kernel  of  sound  pressure  square  with  respect  to  sound  speed  perturbation  between  the  source  and  a 
receiver  underneath  the  second  wave  crest,  and  (c)  the  sound  pressure  sensitivity  and  its  dynamics 
following  the  internal  wave  motion.  The  sensitivity  kernel  has  been  converted  to  be  per  wave  height 

through  a  chain  rule.] 


7 


e  w  (s/w)  ed'< 


—  Sound  Dressure  field 


Distance  (kin) 


04  Nov  2009  09:00:00 
Depth-integrated  Sound  Intenisty 


c  (m/s) 

1520 

1510 

39  05 

1500 
H90 

39 

1480 


40 

-60 

-80 

-100 


38  95 


38.9 


Longitude 


dB 


-10 


-20 


-30 


40 


50 

-60 


-70 


-80 


Figure  2:  Integrated  numerical  simulations  of  ocean  dynamics  and  acoustics  on  the  slope  of  the 

New  Jersey  Shelf 

[(a)The  bathymetry ,  (b)  the  original  down-slope  propagation  pattern,  and  (c)  change  of propagation 
pattern  to  water-borne  ducting  due  to  water  column  variability.  The  beam  pattern  on  the  depth 
integrated  sound  intensity  plots  on  the  right  shows  the  3-D  focusing  of  sound  by  the  small  canyons 

on  the  slope.] 
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Figure  3:  Preliminary  acoustic  data  and  model  analyses  in  the  engineering  test  conducted  on  the 

New  England  Mud  Patch. 

[(a)  and  (b)  Sound  propagation  models  in  two  different  sediment  types  of  environments.  The  sound 
source  in  this  model  study  is  placed  close  to  the  seafloor.  In  the  panel  (a),  the  sediment  layer  is 
modeled  as  fine-granted  sediment  (mud)  with  compressional  wave  speed  less  than  the  sound  speed 
in  the  bottom  of  water  column.  In  the  panel  (b),  the  sediment  layer  is  modeled  as  sand  with  higher 
compressional  wave  speed.  One  can  see  that  the  acoustic  energy  can  be  trapped  in  the  lower  speed 
mud  layer,  (c)  and  (b)  Preliminary  broadband  pulse  data  and  model  comparison.  In  the  panel  (c), 
the  spectrograms  of  received  pulses  on  four  hydrophones  close  to  the  bottom  are  shown,  along  with 
time  series  plots  in  the  bottom  panel.  The  broadband  source  was  made  available  by  D.  Knobles  of 
ARL:UT,  and  the  hydrophone  array  data  was  provided  by  the  ocean  acoustics  group  at  WHOI.  In 
the  panel  (d),  results  from  a  broadband  sound  propagation  model  are  shown.  The  sediment  layer  in 
the  model  is  made  to  be  mud.  Despite  the  very  low  frequency  ground  waves,  this  preliminary  model 
adequately  reproduces  the  pulse  dispersion  over  a  broad  frequency  range.] 
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Figure  4.  Ambient  noise  spectrum  level  during  container  ships  passing  by  the  acoustic  engineering 

test  site  on  the  New  England  Mud  Patch. 

[The  hydrophone  deployment  was  5.5  days,  and  the  noise  spectrum  level  curves  within  a  2  hour  time 
window  around  every  closest  approach  point  of  the  container  to  the  hydrophone  array  are  shown. 
Total  35  containers  were  detected  acoustically,  and  the  total  analyzed  time  for  the  ambient  noise  was 
65  hours  long.  The  black  lines  denote  the  mean  and  one  standard  deviation  curves.  The  color 
curves  are  the  average  levels  of  heavy,  moderate  and  light  shipping  noise  in  deep  ocean  (11).] 
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